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Abstract 

The photoproduction of L)*±(2010) mesons associated with a leading neutron 
has been observed with the ZEUS detector in ep colhsions at HERA using an 
integrated luminosity of 80 pb~^. The neutron carries a large fraction, xl > 0.2, 
of the incoming proton beam energy and is detected at very small production 
angles, 6n < 0.8 mrad, an indication of peripheral scattering. The D* meson 
is centrally produced with pseudorapidity |?7| < 1.5, and has a transverse mo- 
mentum pt > 1.9 GeV, which is large compared to the average transverse mo- 
mentum of the neutron of 0.22 GeV. The ratio of neutron-tagged to inclusive 
D* production is 8.85 ± 0.93(stat.)l°;^^(syst.) % in the photon-proton center-of- 
mass energy range 130 < W < 280 GeV. The data suggest that the presence of 
a hard scale enhances the fraction of events with a leading neutron in the final 
state. 
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1 Introduction 



Events containing a leading neutron have been studied in ep collisions at HERA [1-4]. 
The neutrons carry a large fraction of the incoming proton beam energy, xl > 0.2, and 
are produced at very small scattering angles, 9n < 0.8 mrad, indicative of a peripheral 
process. 

The small transverse momenta (pr) which characterize leading baryon production pro- 
cesses imply a soft scale, which means that a non-perturbative approach is required to 
model such events. Particle-exchange models within Regge theory [5], in particular the 
one-pion-exchange model (OPE) [6-8], are often applied to describe leading neutron pro- 
duction. Charm production, in contrast, can be used to investigate parton dynamics 
because the charm-quark mass provides the hard scale necessary to ensure the applica- 
bility of perturbative Quantum Chromodynamics (pQCD). Therefore the study of charm 
production in events with a leading neutron gives information on the interplay between 
soft and hard scales. 

This letter presents measurements of D*"^ photoproduction associated with a leading neu- 
tron. Differential cross sections and ratios to inclusive D*^ photoproduction are reported. 
These results extend previous ZEUS studies [1,3,4] of leading-neutron production in dijet 
and inclusive photoproduction and deep inelastic scattering. 

2 Charm and neutron production 

An important process in charm photoproduction is Boson-Gluon Fusion (BGF). At leading 
order this corresponds to the direct component where the photon couples directly to a 
high-transverse-momentum cc pair which interacts with a gluon from the proton. Another 
contribution to the cross section comes from the resolved component, where the photon 
acts as a source of partons, interacting with the proton mostly via charm excitation 
processes, cq cq and eg — > eg [9]. 

The mechanism for leading neutron production is not well understood. In the following 
subsections some models for neutron production are briefly discussed. 

2.1 Fragmentation models 

In fragmentation models of partons into hadrons, such as the cluster model [10] or the 
Lund string model [11], a certain fraction of neutrons is expected in the final state. In 
this case the leading neutrons are produced by fragmentation of the proton remnant using 
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the same mechanism as is used for the other final state hadrons. Such models predict a 
softer xl distribution than that measured [1-4]. 



2.2 One-pion-exchange model 

Previous studies have shown that particle-exchange models [6-8,12-15] describe data on 
leading neutron production both at HERA [1-4] and at hadroproduction experiments 
[16-25]. In such models the transition amplitude for p — > n is dominated by OPE and the 
electroproduction cross section can be written as the convolution of a function describing 
the splitting of a proton into a nn system, i.e. the pion fiux factor fn/pixi, t), and the evr 
cross section: 

*K-M. (1) 

where t is the squared four-momentum transfer at the proton vertex, s' = s{l —xl) is the 
squared center-of-mass energy of the en system and s is that of the ep. 

The fiux factors found in the literature can be expressed in general as [8, 15] 

where gp^p/^in) ~ 14.5 is the pup coupling constant, m,r is the pion mass and a(t) is 
defined below. The form-factor F{xL,t) accounts for the finite size of the nucleon and 
pion. Examples of fiux factor parametrizations are: 



F{xL,t) = exp 



a{t) = 0, (3) 



/i [26]: 

(1-^l)J 

where R = 0.6 GeV~^ and F{xL,t) is the light-cone form factor. 
/2 [8]: 

F{xL,t) = 1, ait)=a^{t) (4) 
where aT,{t) ~ t (with t in GeV^) is the Regge trajectory of the pion. 
/3 [27]: 

F{xl, t) = exp[6(t - niDl a{t) = a^{t) (5) 
where h = 0.3 GeV~^ and F{xL,t) is the exponential form factor ; 
h [28]: 

(I _ ^2 /X2\ 

^(^^^^)= \l_//^./ ^ "W = (6) 
where A = 0.25 GeV and F{xL,t) is the monopole form factor. 
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The term a'^'^ in Eq. ^ involve the parton distribution in the pion. Charm production in 
association with a leading neutron is potentially sensitive to the gluon content of the pion 
in OPE models via the BGF process. Parametrizations for the pionic parton distribution 
function (PDF) available in the literature were obtained by performing fits to t^N scatter- 
ing data, assuming some parametrization for the nucleon structure function. Examples of 
such parametrizations are those by Owens [29] which come from fits on J/\& and dimuon 
production data. The more recent GRV parametrizations [30] assume a valence-like struc- 
ture for the pion at a certain low scale. This distribution is dynamically evolved and the 
results combined with the constraints imposed by prompt photon production data on the 
pionic gluon density. 



2.3 Rescattering effects 

From the factorization hypothesis expressed in Eq. it is expected that the ratios, r, 
of neutron-tagged to inclusive cross sections for different electroproduction processes are 
about the same. Most of the dependence of the cross sections on the kinematics of the 
processes cancels; remaining differences can be attributed to differences between the pion 
and proton energies and their PDFs. However, larger differences than these may arise 
from neutron absorption, which can occur through rescattering of the neutron on the 
exchanged photon [14,31]. With increasing size of the virtual photon more rescattering 
may be expected. 

Inclusive phot opro duct ion cross sections are well described by vector meson dominance 
models [32-34], where the dipole associated with the photon is of hadronic size. In 
dijet photoproduction, the presence of the hard scale given by the transverse energy of 
the jets implies smaller dipole size. In the infinite- momentum frame the smaller dipole 
size corresponds to the enhancement of the direct photon component at high transverse 
energy Et- Such an enhancement has been observed at HERA [35]. Absorptive effects 
in dijet photoproduction, therefore, are expected to be smaller than those in inclusive 
photoproduction. In the case of charm photoproduction, an additional hard scale is 
provided by the mass of the charm quark. Previous ZEUS measurements have shown that 
requiring the presence of charm further suppresses the resolved component [36] compared 
to inclusive dijet photoproduction. Therefore rescattering in D* photoproduction may be 
further suppressed in comparison to dijet photoproduction. In deep inelastic scattering 
(DIS), at sufficiently high photon virtuality, little rescattering should occur. Within such 
a picture, therefore, the ratios are expected to have the following relationship: 

j,7P ^jj ^D* ^ ^DIS 
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3 Experimental conditions 



The integrated luminosity of 80.2 ± 1.8 pb used for this measurement was collected at 
the ep collider HERA with the ZEUS detector during 1998 - 2000, when HERA collided 
27.5 GeV electrons or positrons^^ with 920 GeV protons, giving a center-of-mass energy of 
318 GeV. 

A detailed description of the ZEUS detector can be found elsewhere [37]. A brief outline of 
the components that are most relevant for this analysis is given below. Charged particles 
are tracked in the central tracking detector (CTD) [38], which operates in a magnetic 
field of 1.43 T provided by a thin superconducting solenoid. The central tracking detector 
consists of 72 cylindrical drift chamber layers, organized in nine superlayers covering the 
polar-angle^ region 15° < 6* < 164°. The transverse-momentum resolution for full-length 
tracks is a{pT)/pT = 0.0058pT © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (GAL) [39] consists of three parts: 
the forward, the barrel and the rear calorimeters. Each part is subdivided transversely into 
towers and longitudinally into one electromagnetic section and either one (in rear) or two 
(in barrel and forward) hadronic sections. The smallest subdivision of the calorimeter is 
called a cell. The calorimeter energy resolutions, as measured under test-beam conditions, 
are a{E)/E = 0.18/a/E' for electrons and a{E)/E = 0.35/a/E for hadrons, with E in 
GeV. 

The forward neutron calorimeter (FNG) [40,41] was installed in the HERA tunnel at 6* = 
degrees and at Z = 106 m from the interaction point in the proton-beam direction. It is a 
lead-scintillator calorimeter which is segmented longitudinally into a front section, seven 
interaction lengths deep, and a rear section, three interaction lengths deep. The front 
section is divided vertically into 14 towers, allowing the separation of electromagnetic and 
hadronic showers from the energy-weighted vertical width of the showers. The energy 
resolution for neutrons, as measured in a beam test, is a{En)/En = 0.65/a/E^, with neu- 
tron energy. En, in GeV. Three planes of veto counters are used to reject events in which 
particles had interacted with the inactive material in front of the FNG. Magnet apertures 
limit the FNG acceptance to neutrons with production angles less than 0.8 mrad, which 
corresponds to transverse momenta pr < EnO^^^^ = 0.74xl GeV. The mean value of pr 
for the data is 0.22 GeV. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > e'-fp, 
^ Hereafter, both e'^ and e~ are referred to as electrons. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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where the photon was measured with a lead-scintillator calorimeter [42,43] located at 
Z = -107 m. 

4 Kinematics 

The kinematics of photoproduction at HERA are specified by the photon virtuality, Q^, 
and the photon-proton center-of-mass energy, W. The electron-proton center-of-mass 
energy, ^/s, is related to W by W"^ = ys where y is the fraction of the electron beam 
energy carried by the photon in the proton rest frame. 

To describe the process ep — > e'D*^nX, four additional variables are used: two for the 
neutron and two for the charmed meson. They are: 

• {xL,On), the fractional energy and production angle of the produced neutron; only 
about half of the data have a On measurement, therefore all results discussed here are 
integrated over this variable up to the maximum accepted angle of 0.8 mrad; 

• {PT,v)y the transverse momentum and pseudorapidity of the produced D*^ meson. 

The measurement was performed in the following kinematic region: < 1 GeV^, 
130 <W < 280 GeV, \v{D*)\ < 1.5, Pt{D*) > 1.9 GeV, xl > 0.2 and ^„ < 0.8 mrad. 

5 Event Selection 

5.1 Trigger 

A three-level trigger system was used to select events online [37,44]. The selection was 
based on energy deposits, tracking and event timing. The FNC was not used in the 
trigger. 

5.2 Photoproduction selection 

Photoproduction events were selected offline using cuts based on the reconstructed pri- 
mary vertex position, CAL energy deposits and the reconstructed tracks of charged parti- 
cles. Events with a well-identified electron candidate in the CAL were removed. It was re- 
quired that '^i{Ei—pz,i) > 7 GeV, where the sum runs over all CAL cells and pz,i is the Z 
component of the momentum vector assigned to each cell of energy Ei. Tracking and CAL 
information was combined to form energy fiow objects (EFOs) [45,46]. A cut was made on 
the Jacquet-Blondel [47] estimator of W"^, Wj^ = i/jbs, where yjB = J^ii^i ~ Ez,i)/2Ee, 
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and Ez^i = EiCosOi] Ei is the energy of EFO i with polar angle 6i with respect to the 
measured Z- vertex of the event. The sum runs over all EFOs. It was required that 
WjB < 265 GeV. These cuts correspond to a true W range of 130 < W < 280 GeV and 
g2 < 1 GeV^ with the median ^ IQ-'^ GeV^. 

5.3 D*^{2010) reconstruction 

The inclusive charm sample was selected by identifying events containing a charmed me- 
son. The D*^ selection cuts are based on the decay channel: D*^ {D^ K~7r~^)7r'^ 
(+ charge conjugate), where vr,; indicates the "slow" pion [48]. Only tracks assigned to 
the primary event vertex and with hits in at least three superlayers of the CTD were 
considered. The combinatorial background was reduced and the kinematic phase space 
defined by requiring: the transverse momentum of the kaon and pion candidates to satisfy 
Pt{K) > 0.45 GeV, Pt{t^) > 0.45 GeV and Pri'^s) > 0.12 GeV; the transverse momen- 
tum of the D*^ to be greater than 1.9 GeV and the pseudorapidity of the D*^ to satisfy 
\viD*^)\ < 1-5. 

Since no particle identification was performed, the K and vr masses were alternately at- 
tributed to the decay products of the candidate D° meson. Those candidates that 
had an invariant mass between 1.80 and 1.92 GeV were required to have a mass differ- 
ence AM = M{K7i7is) - M{Kti) between 0.1435 and 0.1475 GeV. The combinatorial 
background was estimated from the mass-difference distribution for wrong-charge combi- 
nations, in which both tracks forming the candidates have the same charge and the 
third track has the opposite charge. 

5.4 Neutron reconstruction 

Events with a leading neutron were selected from the inclusive charm sample by requiring 
a large energy deposit [En > 184 GeV) in the ENG. Protons, photons, and neutrons are 
separated by their position in the detector, as well as by the shower width. Scattered 
protons are deflected by the HERA magnets and strike the top part of the ENG. Pho- 
tons can be identified and removed from the sample because the transverse spread of 
electromagnetic showers is much less than that of hadronic showers. 

Events with particles that started to shower before reaching the ENG were removed by re- 
quiring that the scintillator veto counter had an energy deposit below that of a minimum- 
ionizing particle. Events with wide showers, inconsistent with originating from a single 
high-energy hadron, were removed. 
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5.5 Final event sample 



The AM distribution for the neutron-tagged sample is shown in Fig. ^ along with the 
wrong-charge combinations. A prominent D*"^ signal is observed. The signal observed in 
the M{Ktt) distribution for events within the mass window 0.1435 < AM < 0.1475 GeV 
is shown as an inset. 

After the wrong-charge background subtraction, 298 ± 31 mesons were found. 

The same background subtraction procedure, applied to the inclusive D*"^ sample, gave 
14 743 ± 253 events. 

6 Monte Carlo simulation and acceptance corrections 

A GEANT-based [49] Monte Carlo (MC) simulation was used to calculate selection efficien- 
cies and correction factors for the charmed meson. Three different event generators were 
used: Rapgap 2.08/06 [50] for evaluating the nominal corrections, Herwig 6.301 [51] and 
Pythia 6.156 [52] as systematic checks. Rapgap and Pythia use the Lund string model 
for hadronization. Herwig uses a cluster model. The events generated with Rapgap for 
acceptance calculations were produced using OPE for the production of the leading neu- 
tron, with the pion flux factor from Eqs. (j2I) and Q and the GRV parametrization [30] for 
the pion PDF. Inclusive Rapgap, employing the Lund string model instead of OPE for 
neutron production, was produced for comparisons to the final measurements. The leading 
neutron is also produced via the Lund string model in Pythia. In Herwig, it is pro- 
duced via the cluster model. The proton PDFs parametrizations used were CTEQ5L [53] 
for Pythia and Herwig, and CTEQ4D [54] for inclusive Rapgap. The photon PDF 
GRV-G LO [55] was used in Pythia and Herwig, and GRS LO [56] in Rapgap. The 
mass of the charm quark was set to 1.5 GeV. The fraction of c quarks hadronizing to a D* 
meson was set to /(c D*) = 0.235 [57]. Both direct and resolved photon processes for 
charm production were generated, in proportion to their predicted cross sections. Effects 
of neutron rescattering were not taken into account in the simulation. 

For all the MC samples used to evaluate the acceptances, events with at least one D*^ 
decaying in the appropriate decay channel were selected and passed through the ZEUS 
detector and trigger simulations as well as the event-reconstruction package. 

Since the D* and neutron were independently detected and their kinematics largely un- 
correlated, the acceptances for the two particles factorize. The selection efficiencies and 
correction factors for the neutron calculated for previous analyses [4] were used. The 
overall acceptance of the FNC, which includes the beam-line geometry and the angular 
distribution of the neutrons, is about 25% for neutrons with xl > 0.2 and On < 0.8 mrad. 
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The differential cross section for D*^ photoproduction associated with a leading neutron 
was evaluated in terms of a given variable Y as da/dY = N/ (Afnc ■ Ao* -B- AF) where N 
is the number of D* found in the final sample in a bin of size AF, A-p^Q is the acceptance 
for the neutron detection in the FNC, Ad* is the acceptance for the D* reconstruction 
and B is the branching ratio for the selected decay mode. A value ol B = 2.57% [58] was 
used. 

7 Systematic uncertainties 

For the -D*^ measurement the major sources of systematic uncertainty are listed below 
(the relative uncertainty on D*"^ acceptance is shown in parentheses): 

• the selection of photoproduction events and D*^ candidates. Variations were made in 
the PFjB (+4.5 %) and Z-vertex (+JJ %) cuts; 

• the Pt of the pion and kaon cuts. These were varied according to their resolutions 

(ti:? %); 

• variation of the mass windows. The AM window used for the extraction of the D*"^ 
was widened symmetrically by 0.5 MeV. The M{D^) window was widened and reduced 
symmetrically by 5 MeV {t?^^ %); 

• the AM region for the normalization of the wrong-charged combinations. This was 
changed from 0.15 - 0.165 GeV to 0.15 - 0.163 GeV (+2.4 %); 

• the MC model dependence. Herwig (-4.8 %) and Pythia (+2 %) were used instead 
of Rapgap; 

• the fraction of resolved photon events in the MC was lowered by 20% and raised by 
10 % (1?| %); 

• the CAL energy scale. This was varied within its uncertainty of ±3% (±0.9 %). 

An extensive discussion of the systematic effects related to the neutron measurement is 
given elsewhere [4]. Here, the major sources of systematic uncertainty and their effect on 
the FNC acceptance (shown in parentheses) is listed: 

• the uncertainty in the angular distribution of the neutrons (±4% for xl < 0.82, 
±7% for XL > 0.82); 

• the uncertainty in the overall FNC energy scale of ±2% (less than 4% effect for 
XL < 0.82, tlt% for XL > 0.82 ); 

• the normalization uncertainty arising from proton-beam-gas interactions overlapping 
with photoproduction events, the uncertainty in the amount of dead material in the 
beam line, and the uncertainties from the veto cuts (±5%). 
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All above errors were added in quadrature separately for the positive and negative vari- 
ations to determine the overall systematic uncertainty. The overall normalization has 
additional uncertainties of 2.2% due to the luminosity measurement and 2.5% due to the 
knowledge of branching ratios. 

Sources of systematic uncertainty in the ratio measurement were studied in a similar 
manner to those for the cross-section measurements. There is a cancellation between 
the common systematic uncertainties originating from the selection of inclusive photo- 
production events, the selection of D*^ candidates and the background estimation. The 
remaining contributions are those from the model dependence of the acceptance correc- 
tions used in the evaluation of the inclusive Z)*^ photoproduction cross sections and from 
the neutron measurement uncertainties. 

8 Results 

The integrated cross section for the reaction ep e'D*^nX in the kinematic region 
< 1 GeV^, 130 <W < 280 GeV, \v{D*)\ < 1.5, Pt{D*) > 1.9 GeV, xl > 0.2 and 
On < 0.8 mrad is 

2.08 ± 0.22(stat.)+°;j^(syst.) ± 0.05(B.R.) nb, 

where the final uncertainty arises from the uncertainty of the branching ratios for the 
D* and D^. The luminosity uncertainty was included in the systematic uncertainty. The 
predicted cross sections from the models are 3.0 nb for Herwig, 4.6 nb for Pythia, 2.6 nb 
for inclusive Rapgap and 2.0 nb for Rapgap with OPE. 

Table 1 and Fig. |21 show the differential cross sections for neutron-tagged D*^ production 
as a function of W, pt{D*) and ri{D*). The differential cross section as a function of 
is shown in Table 1 and Fig. El 

The inclusive D*^ cross section was measured and found to agree with previous measure- 
ments [36] in a similar kinematic range. Table |21 and Fig. |3] show the measured ratios of 
neutron-tagged to inclusive D*^ production as a function of different kinematic variables. 
Over the whole measured kinematic range the ratio is 

r^* = 8.85 ± 0.93(stat.)+|]:^^(syst.) % 

which is shown superposed on Fig. EJ The P^r degree of freedom with respect to the 
overall ratio are 0.27, 1.65 and 0.09 for the W, pt{D*) and rj^D*) distributions. Within 
the experimental uncertainties neutron-tagged D*^ production is compatible with being 
a constant fraction of inclusive D*^ production, independent of the D*^ kinematics. 



9 



9 Discussion 



The experimental results in Figs. El and Ok are compared to the predictions, normalized to 
the data, of the MC models Rapgap with OPE for leading-neutron production, Herwig, 
Pythia and inclusive Rapgap. In Fig. |2l the agreement in shape between data and MC 
models is satisfactory in all cases. Fig. |2K; however, shows that only Rapgap with OPE 
agrees with the measured neutron energy distribution seen in the data. A similar result 
was obtained in the study of neutron-tagged dijet photoproduction [3]. 

In principle the data allow the pion PDF to be probed in the range of the parton frac- 
tional momenta 10"'^ < < 10"^. Figure [Hb shows the differential cross section in xl 
compared to the predicted cross sections from Rapgap with OPE, based on different 
parametrizations for the pion structure function: GRV set 1 and Owens sets 1 and 2. The 
data have little sensitivity, as in the case of neutron-tagged dijet photoproduction [3], to 
the choice of the pion structure function. Even with an extreme choice of a pion structure 
function, e.g. a completely fiat gluon distribution, or a parametrization identical to the 
proton structure function, little variation is seen in the predictions. 

Figure |St shows the data compared to the predictions of Rapgap with OPE for the four 
flux factors discussed in Section 12.21 All Rapgap distributions are normalized to the 
data and the resulting normalization factors are given in the figure. Fluxes /i and /s give 
similar results, being compatible with the data both in shape and normalization. The 
fluxes /2 and are disfavored by the shape of the data, and the latter predicts cross 
sections almost three times smaller than data. 

The independence of r^* on W, pt{D*) and ri[D*), shown in Fig. EJ supports the hy- 
pothesis of vertex factorization, in agreement with previous studies of leading-neutron 
production [1,3,4]. The predicted ratios from the models are: 0.24 for Herwig, 0.29 for 
Pythia, and 0.18 for Rapgap. 

The ratio for charm production, r^* , is in agreement with the analogous ratio previously 
measured for neutron-tagged DIS, > 4 GeV^ [4], 

r°^s = 8.0 ± 0.5 %, 

but lies above the ratio previously measured for neutron-tagged inclusive photoproduction 
atW = 207 GeV [4], 

r^P = 5.7 ± 0.4 %, 

as expected from the rescattering effects within the OPE model (see Section f2.3j) . The 
errors quoted are the quadratic sum of the statistical and systematic uncertainties. 

For photoproduced dijets the neutron-tagged to inclusive ratio has only been measured 
for xl > 0.49. In this kinematic region, the measured D*^ ratio and their corresponding 
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ratios previously measured for DIS, inclusive photoproduction [4], and photoproduction 
of dijets with transverse energy i?;^* > 6 GeV [3] are: 



r^*(xL > 0.49) = 6.55±0.76(stat.)^o;45(syst.) % 

r°^^(xL > 0.49) = 5.8 ±0.3% 

rjj(xL>0.49) = 4.9 ±0.4% 

r^P(xL > 0.49) = 4.3 ± 0.3 %. 

The results are compatible with the rescattering hypothesis described in Section 2.3. 



10 Summary 

The photoproduction of D*"^ mesons associated with a leading neutron has been studied 
in ep interactions at HERA in the kinematic region <1 GeV^, 130 < < 280 GeV, 
\ri{D*)\ < 1.5, Pt{D*) > 1.9 GeV, 6^ < 0.8 mrad and xl > 0.2. The Monte Carlo mod- 
els Rapgap, Herwig and Pythia give a satisfactory description of the D*^ kine- 
matics, but only Rapgap with one-pion exchange satisfactorily describes the leading- 
neutron energy distribution. The results show sensitivity to the choice of pion flux factor. 
The ratio of neutron-tagged D*^ photoproduction to inclusive D*^ photoproduction is 
r^* = 8.85 ± 0.93(stat.)^|]:^?(syst.)%. 

The ratio of neutron-tagged D*^ photoproduction to inclusive D*^ photoproduction is 
constant as a function of W, pt{D*) and rj{D*), in agreement with the hypothesis of 
vertex factorization. This ratio is consistent with the analogous ratio in deep inelastic 
scattering, but both are about 30% higher than the corresponding ratio for inclusive 
photoproduction, suggesting that the presence of a hard scale enhances the fraction of 
events with a leading neutron in the final state. 
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W Range ( GeV) 


da/dW ± (stat.) ± (syst.) (nb/GeV) 


130 - 160 


0226 ± 0054'^nnnol 

KJ • KJ KJ — 1 — KJ • KJ KJ i- r\ flMQ 


1 60 - 1 88 

X \J\J ± KjKj 


0214 + 0036+S-2S^I 

KJ .KJ jLi Xt: _1_ V 7 . V 7 V7 1. J V J rj nnfl^. 


188 - 226 


0111 ± 0026+nnf!n? 

KJ * KJ -L -L -L — 1 — KJ t KJ KJ KJ f\ lllll}/ 


226 - 280 


0.0080 ± 0.0017j:ooo?7 


Vt(D*^) Ranee f GeV) 


da/dvT(D*^) ± fstat ) ± fsvst ) fnb/GeV) 


1.9 - 2.3 


1.991 ± 0. 531^0 421 


2.3 - 2.73 


1 089 ± 261'*"n}fQ 

_l_ • KJ K^ \J — 1 — KJ t KJ -L I 1 1 


2.73 - 3.8 


364 ± 090"'"nm c 

KJ t '-.J KJ -L — 1 — KJ t KJ \J KJ / 1 1 J 1 V 


3.8 - 15 


0.038 ± 0.004+0 004 


n( D*^) Ranee 


doldn{D*^\ ± fstat ) ± fsvst ) fnb) 




914 ± 126"'"oS?I 


f— 72") - f— 15') 


858 + 147+1^-irr 

u.OtJO in u. I'll _o.055 


(-0.15) - (+0.42) 


0.665 ± 0.161+nnR^ 

V / . vy w -v_/ _i_ vy . -1- vy -L _|j 


(+0.42) - (+1.5 ) 


0.436 + 0.132+°:?^^ 


xl Range 


da/dxL ± (stat.) + (syst.) (nb) 


0.2 - 0.46 


2.18 ± 0.48l|^:2o 


0.46 - 0.64 


3.70 ± 0.64l°J^ 


0.64 - 0.82 


4.29 ± 0.65lSi^ 


0.82 - 1.0 


0.456 + 0.381^|];J^^^ 



Table 1: Values of the differential cross sections for neutron-tagged D photo- 
production (Q^ < 1 GeV^ andOn < 0.8 mrad) with respect to W , pt{D*^), rj{D*^) 
and xl- 
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W Range ( GeV) 


r^' ± (stat.) ± (syst.) 


130 - 160 
160 - 188 
188 - 226 
226 - 280 


0.089 ± 0.0221°:°°^ 

0.101 ±o.oi7l|J:[{^^ 

0.075 ± ^mStomi 
0.090 ± 0.019t|];[J^^ 


Pt{D*^) Range (GeV) 


r^' ± (stat.) ± (syst.) 


1.9 - 2.3 
2.3 - 2.73 
2.73 - 3.8 
3.8 - 15 


0.137 ±0.038l|]:[{|^? 
0.105 ±0.026l|J:[{^^ 
0.058 ± 0.014t°:°J^ 
0.087 + 0.010+°;°°^ 


■r]{D*^) Range 


r^' ± (stat.) ± (syst.) 


(-1.5 ) - (-0.72) 
(-0.72) - (-0.15) 
(-0.15) - (+0.42) 
(+0.42) - (+1.5 ) 


0.093 ± 0.013^|];[^^^ 

0.097 ± o.oi6l|J:[{^I 

0.085 ± 0.0211°:°°^ 
0.080 ± 0.024+°;°°^ 



Table 2: Values of the ratios of the differential cross sections for neutron-tagged to 
inclusive D*^ photoproduction (Q"^ < 1 GeV^, xl > 0.2 and 6n < 0.8 mrad) with 
respect to W , pt{D*^) and ■r]{D*^). 



16 



ZEUS 

• ZEUS 1998-2000 
Xl > 0.2, < 0.8 mrad 
Wrong-charge background 
N(D* ^ Knn) = 298±31 




0.145 



0.15 0.155 0.16 0.165 

M(Knn) - M(K7i) (GeV) 



Figure 1: The data points show the neutron-tagged AM distribution for right- 
charge track combinations. The solid line shows the wrong-charge combinations 
normalized to the right-charge combinations in the region 0.15 < AM < 0.165 GeV 
outside the D*^ mass window 0.1435 < AM < 0.1475 GeV which is shown shaded. 
The signal observed in the M{Kti) distribution for events within the mass window 
0.1435 < AM < 0.1475 GeV is shown as an inset. The solid line also shows the 
wrong-charge combinations normalized to the right-charge combinations as before. 
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Figure 2: The points show the differential cross sections for neutron-tagged D*^ 
production as a function ofW, p^iD*) andrj^D*) for xl > 0.2 and 6^ < 0.8 mrad. 
The error bars displayed on the plots denote the statistical uncertainty (inner) and 
the quadratic sum of the statistical and the systematic uncertainties (outer). The 
uncertainties due to the luminosity measurement and the branching ratios are not 
shown. The predictions of Monte Carlo models normalized to the data are also 
shown. 
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Figure 3: The points show the differential cross sections for neutron-tagged D*^ 
production as a function ofxi. The histograms show the predictions of Monte Carlo 
models (a) Rapgap with OPE (solid histogram), Herwig (dashed), Pythia (dot- 
ted), and inclusive Rapgap (dashed- dotted); (b) Rapgap with OPE and pion PDF 
parametrizations from GRV set 1 (solid), Owens sets 1 (dashed) and 2 (dotted); 
(c) Rapgap with OPE and flux factors /1-/4 from Eqs. The error bars 

displayed on the plots denote the statistical uncertainty (inner) and the quadratic 
sum of the statistical and the systematic uncertainties (outer). The uncertainties 
due to the luminosity measurement and the branching ratios are not shown. All 
distributions are normalized to the data. The numbers in (c) are the normalization 
factors. 
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• ZEUS 1998-2000, 9^ < 0.8 mrad, 0.2 




140 160 180 200 220 240 260 280 

W (GeV) 




Figure 4: The ratio of neutron-tagged D*^ production to inclusive D*^ production 
as a function of W , p^{D*) and t]{D*) for xl > 0.2 and On < 0.8 mrad. The 
error bars displayed on the plots denote the statistical uncertainty (inner) and the 
quadratic sum of the statistical and the systematic uncertainties (outer). The line 
superposed on the figures shows the overall ratio of neutron-tagged D*^ to inclusive 
D*^ events. 
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